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ABSTRACT 


In  understanding  the  complexity  of  the  fracture  mechanics  of  brittle  matrix 
composites,  knowing  reliable  values  of  the  mechanical  properties  of  the  fiber-matrix 
interface  are  important.  Since  all  experiments  designed  to  evaluate  the  mechanical 
properties  of  the  fiber-matrix  interface  rely  on  curve  fitting  experimental 
measurements  to  analytical  results,  accurate  physical  modeling  of  these  tests  is 
imperative. 

In  this  final  report,  two  common  tests,  namely  the  push-in  test  and  the  slice 
compression  test,  were  analytically  simulated.  In  the  push-in  test,  the  effect  of 
extrinsic  factors  such  as  indentor  shape  and  size,  and  push-through  radius  hole,  and 
intrinsic  factors  such  as  transverse  isotropy  of  fibers  on  load-displacement  curves 
was  studied.  All  factors  except  the  push  through  hole-radius  were  found  to  affect 
the  load-displacement  curves. 

In  the  slice  compression  test,  comparison  was  made  for  the  values  of 
maximum  protrusion,  residual  protrusion  and  debond  lengths  for  a  simple  shear- 
lag  analysis  and  a  finite  element  analysis  model.  Large  differences  ranging  from 
15%  to  70%  were  foimd  by  using  the  two  different  models.  Hence,  the  importance 
of  accurate  modeling  of  the  test  cannot  be  ignored.  It  also  shows  that  there  is  a  need 
for  further  study  in  this  area  to  synergize  the  experimental  and  analytical  tools 
required  to  find  the  mechanical  properties  of  the  fiber-matrix  interface.  The  study 
also  gives  tools  to  an  experimentalist  for  designing  reliable  experiments  for 
evaluating  reliable  values  of  the  mechanical  properties  of  the  fiber-matrix  interface. 
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INTRODUCTION 


The  advantages  of  using  ceramics  are  their  high  refractoriness  and  strength, 
low  density,  excellent  corrosion  and  oxidation  resistance,  and  low  cost.  However, 
their  brittle  nature  and  low  fracture  toughness  makes  them  prone  to  catastrophic 
failure.  A  promising  method  to  overcome  these  shortcomings  is  by  reinforcing 
ceramics  by  aligned  continuous  fibers,  such  as,  silicon  carbide  and  carbon.  The 
resulting  composite  materials  are  called  ceramic  matrix  composites  (CMCs),  and 
have  increased  fracture  toughness  and  fail  more  gracefully.  However,  they  are 
highly  anisotropic  and  fail  in  a  complex  manner,  and  often  do  not  follow  the  basic 
concepts  of  fracture  mechanics.  So,  the  basic  understanding  and  modeling  of  the 
mechanics  of  these  materials  are  essential  for  the  development  of  these  materials. 

In  this  study,  two  common  tests  used  to  find  the  values  of  the  mechanical 
properties  of  the  fiber-matrix  interface  were  theoretically  simulated.  Accurate 
modeling  of  these  tests  is  especially  critical  as  experimental  measurements  are  curve 
fitted  to  analytical  models  to  evaluate  ihe  mechanical  properties  of  the  fiber-matrix 
interface.  For  example,  in  one  of  the  two  tests  simulated  here,  namely,  the  push-in 
test,  the  experiment  measures  the  load  vs.  displacement  data.  Simplified  analytical 
models  develop  load-displacement  curves  as  a  function  of  mechanical  properties  of 
the  constituents,  fiber  volume  fraction,  and  fiber-matrix  properties  -  coefficient  of 
friction  and  residual  stresses.  Knowing  the  mechanical  properties  of  the 
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constituents  and  the  fiber  volume  fraction,  one  can  regress  to  find  the  coefficient  of 
friction  and  residual  stresses  at  the  fiber-matrix  interface.  Although  simple 
analytical  models  are  easy  to  use,  are  they  accurate  for  finding  values  for  the  two 
interface  mechanical  properties?  This  is  the  central  theme  of  this  work  and  is 
discussed  under  Task  1  and  Task  2  of  this  final  technical  report. 

All  of  the  work  done  in  this  grant  is  available  in  thesis  and  refereed  journals. 
The  titles,  abstracts  and  their  references  are  given  in  Appendix  A.  The  cumulative 
list  of  researchers  is  given  in  Appendix  B.  We  give  a  comprehensive  executive 
summary  of  the  significant  work  accomplished  in  this  grant  in  the  next  section. 
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TASKS 


Task  1:  Effect  of  Extrinsic  and  Intrinsic  Factors  on  an  Indentation 

Test 

The  interface  is  a  vital  part  of  the  ceramic  matrix  composite  in  determining 
its  strength  and  toughness.  Hence  quantitative  value  of  the  mechanical  properties 
of  the  interface  between  the  fiber  and  matrix  is  of  great  importance.  However,  all 
experiments  designed  to  find  these  values  rely  on  curve  fitting  the  experimental 
measurements  to  theoretical  models.  Hence  the  reliability  of  the  obtained 
mechanical  properties  of  the  interface  is  dependent  on  the  accuracy  of  the 
theoretical  models. 

One  of  the  tests  conducted  for  the  measurement  of  the  mechanical  properties 
of  the  interface  is  the  push-out  test  (Figure  1).  Several  analytical  models  have  been 
developed  for  the  test  md  they  range  form  simple  analytical  models  based  on  shear- 
lag  assumptions  (Shetty,  1988)  to  more  physically  realistic  finite  element  models 
(Shirazi-Adi,  1992).  However,  the  effect  of  extrinsic  and  intrinsic  factors  on  the 
push-out  test  has  not  been  studied.  These  factors  include  the  following: 

Extrinsic  Factors 

1.  Shape  (smooth,  flat  and  pointed)  of  indentor  (Figure  2) 

2.  Size  (contact  radius  to  fiber  radius  ratio)  of  indentor 

3.  Radius  of  the  push-through  hole 
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Intrinsic  Factors 


1.  Relative  elastic  moduli  of  fiber  and  matrix 

2.  Material  symmetry  (Isotropic  vs.  Transversely  Isotropic)  of  fiber 

Both  analytical  (Figure  3)  and  finite  element  models  (Figure  4)  (Maze,  1983; 
Nike-2D,  1991;  Orion,  1985)  were  used  to  understand  the  effect  of  the  above  factors. 
In  the  analytical  model,  a  semi-infinite  circular  fiber  is  bonded  to  a  half-space 
matrix.  The  fiber  and  matrix  are  assumed  to  have  identical  mechanical  properties 
and  the  interface  either  is  perfectly  bonded  or  follows  the  Coulomb  friction  law. 
The  finite  element  geometry  extends  the  above  model  to  fiber  and  matrix  of  finite 
length  and  radius. 

The  following  conclusions  were  drawn  from  the  various  studies; 

•  The  indentation  test  is  not  affected  by  the  radius  of  the  push-through  hole 

•  The  displacements,  slip  lengths  and  stresses  for  the  flat  indentor  case  are 
more  than  that  in  the  uniform  pressure  case  (Figure  5).  Also,  the  axial 
displacement  between  the  fiber  and  the  matrix  increases  as  the  radius  of  the 
indentor  increases  (Figure  6).  Since  this  is  one  of  the  experimental 
measurements,  the  shape  and  radius  of  indentor  should  be  dealt  in  the  curve 
fitting  of  the  indentation  tests. 

•  Assuming  a  transversely  isotropic  fiber  to  be  isotropic  is  not  suitable  as  the 
load  vs.  fiber-matrix  interface  displacements  vary  substantially  (Figure  7). 
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Task  2:  Simulation  of  the  Slice  Compression  Test 


This  task  also  concerns  the  need  for  accurate  modeling  of  theoretical  models 
for  reliable  calculation  of  mechanical  properties  of  the  fiber-matrix  interface.  In  this 
study,  the  slice  compression  test  (Shafry,  et  al  1989,  Shafry  and  Brandon,  1991) 
(Figure  8)  is  simulated  to  find  the  effect  of  using  approximate  techniques,  such  as, 
shear-lag  analysis  to  model  the  test.  This  is  important  as  slice  compression  test  like 
the  indentation  test  does  not  directly  evaluate  the  mechanical  properties  of  the  fiber- 
matrix  interface.  If  the  analytical  model  used  to  regress  the  experimental  data  to 
find  the  mechanical  properties  of  the  interface  is  questionable,  the  results 
themselves  are  not  reliable  either. 

In  contrast  to  the  push-in  test  of  Task  1,  the  slice  compression  test  can  be  used 
to  extract  interfacial  properties  where  a  large  number  of  fibers  are  tested 
simultaneously.  In  this  test,  a  thin  slice  of  unidirectional  aligned  composite  (about 
2mm  thick)  is  cut  normal  to  the  fiber  direction  and  the  opposite  flat  faces  are 
polished  (Figure  8).  This  slice  is  then  compressed  between  a  soft  top  plate  (typically 
made  of  aluminum  or  copper  alloys)  and  a  hard  bottom  plate  (typically  made  of 
silicon  nitride).  The  loading  on  top  of  the  aluminum  plate  is  slowly  increased  to  a 
peak  value,  and  then  unloaded  completely. 

On  loading  the  specimen,  at  sufficiently  high  loads,  debonding  is  initiated  at 
the  top  surface  of  the  specimen  between  the  fiber  and  matrix  (Figxue  8).  This  causes 
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the  stiffer  component  (usually  the  fibers  in  ceramic  matrix  composites)  to  protrude 
into  the  aluminum  plate.  This  creates  a  negative  impression  on  the  soft  aluminum. 
On  unloading,  the  fibers  relax  partially  back  into  the  matrix,  thus  leaving  a  residual 
protrusion  of  the  fibers. 

Since  aluminiun  is  very  plastic,  the  depth  of  the  indentation  on  the  aluminum 
plate  is  usually  taken  as  equal  to  the  maximum  protrusion  of  the  fibers  at  maximiun 
load.  The  maximum  and  the  residual  protrusion  are  measured  in  this  test  using 
scanning  electron  microscopes. 

The  quantities  thus  obtained  have  to  be  related  to  the  interfacial  parameters, 
namely  the  coefficient  of  friction  and  residual  stresses.  For  this  purpose  a 
mathematical  model  of  the  test  is  required  (Figure  9).  Analytical  studies  have  been 
made  of  the  slice<ompression  test  by  several  authors  (Lu  and  Mai,  1994  and  Hseuh, 
1993).  All  models  are  based  on  shear-lag  analysis.  Hence  there  was  a  need  to 
develop  a  physically  more  realistic  model  such  as  based  on  finite  element  analysis. 

In  this  task,  the  slice  compression  test  loading  is  modeled  in  two  ways:  a)  by 
distributing  the  overall  load  on  the  top  plate  by  using  unequal  but  uniform  loads 
on  the  fiber  and  matrix  as  done  in  shear-lag  models,  b)  by  modeling  the  top  plate 
as  elastoplastic  and  applying  a  uniform  pressure  to  it.  The  maximiun  and  residual 
protrusions,  and  debond  lengths  are  then  compared  between  the  finite  element  and 
shear-lag  analysis  models  for  both  cases  of  loading.  The  effect  of  finite  modulus  of 
the  base  plate  on  the  protrusions  obtained,  and  the  distribution  of  pressure  on  the 
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surface  of  the  fiber  and  matrix  on  loading  the  aluminum  plate  is  also  studied. 

The  conclusions  drawn  from  the  studies  were: 

•  The  three  measurements  -  maximxun  protrusion,  residual  protrusion  and  debond 
length  are  highly  dependent  on  the  modeling  technique.  Differences  occur 
anywhere  from  10  to  70%  between  the  finite  element  and  shear-lag  models  (Figures 
10-12). 

•  Assuming  the  base  plate  as  non-rigid  made  a  difference  of  15%  in  some  of  the  above 
measurements  (Figure  13). 

•  Assuming  uniform  stresses  on  the  top  surface  of  the  fiber  and  matrix  is  not  valid 
(Figure  14).  Even  the  average  stresses  on  the  fiber  and  matrix  found  by  directly 
applying  the  stresses  to  the  top  plate  differed  by  as  much  as  25%  with  the  shear-lag 
models. 

All  the  above  results  show  that  slice  compression  tests  need  to  be  regressed 
to  models  which  are  physically  more  realistic.  Since  the  results  of  the  slice 
compression  test  do  not  directly  give  the  values  of  the  coefficient  of  friction  and 
residual  stresses,  but  by  regressing  it  to  an  analytical  model,  any  trust  in  such  values 
needs  to  be  verified.  Clearly  shear-lag  and  finite  element  models  give  quite 
different  responses.  Although  the  former  model  explicitly  gives  the  values  of  the 
two  mechanical  properties  of  the  interface  from  the  three  experimentally  measured 
values  of  maximum  protrusion,  residual  protrusion  and  debond  length,  the  same 
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cannot  be  said  about  finite  element  models.  So  further  study  is  required  in  this  area 
to  come  up  with  a  reasonably  acceptable  analytical  model  for  regressing  to  find 
these  interface  properties. 

Both  the  above  studies  are  critical  in  the  ongoing  Air  Force  effort  to 
characterize  brittle  matrix  composites  for  use  in  high  temperature  applications  such 
as  jet  engines.  Developing  confidence  in  the  experimental  and  analytical 
characterization  of  brittle  matrix  composites  are  prerequisites  to  a  reliable  and 
optimum  design  of  components  made  of  such  advanced  composite  materials. 
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Title: 

Authors: 

Journal: 

Abstract : 


Effect  of  Extrinsic  and  Intrinsic  Factors  in  Indentation  Tests 
T.  Srinath,  H.  Madanraj,  A.K.  Kaw,  J.  Ye  and  G.H.  Besterfield 
International  Tournal  of  Solids  and  Structures.  Vol.  33,  No.  24,  pp. 
3497-3516,  (1996). 

The  effect  of  extrinsic  and  extrinsic  factors  on  measured  results,  such 
as  load-displacement  curves  and  interfacial  stresses,  from  indentation 
tests  of  composite  materials  is  studied  using  both  analytical  and  finite 
element  models.  The  intrinsic  factors  include  properties  of  the  fiber- 
matrix  interface  and  the  material  symmetry  of  the  fiber  (transversely 
isotropic  or  isotropic).  The  extrinsic  factors  include  the  radius  of  the 
hole  through  which  the  fiber  is  pushed  in,  and  the  size  and  shape  of 
the  indentor.  Out  of  the  above  factors,  only  the  radius  of  the  hole  has 
negligible  effect  on  the  results  of  the  indentation  test. 
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Title: 

Authors: 

Journal: 

Abstract: 


Should  We  Continue  to  Debate  over  the  Change  to  SI 
System? 

A.K.  Kaw  and  M.  Daniels 

Journal  of  Professional  Issues  in  Engineering  Education  and 
Practice.  Vol.  122,  No.  2,  pp.  69-72,  (1996). 

Ninety  five  percent  of  the  world's  population  uses  the  international 
system  of  units  -  Systeme  International  (SI).  However,  we  continue 
to  debate  our  country  converting  to  the  international  system  of  units. 
In  this  paper,  we  have  presented  a  comprehensive  list  of  the 
advantages  of  the  SI  system.  We  then  give  the  reasons  why  we  need 
to  convert  to  the  SI  system  of  imits.  A  few  major  corporations  in  USA 
have  already  made  the  changeover  to  the  SI  system,  but  the  general 
population  is  still  reluctant  to  make  the  change  from  the  US 
customary  imits.  Studying  these  conflicting  issues  simultaneously  is 
important  in  the  implementation  of  the  SI  system,  if  enforced  in  our 
coimtry.  We  have  also  enumerated  the  ways  we  can  progress  toward 
the  changeover.  These  ways  include  learning  from  the  experiences  of 
our  own  failed  effort  of  the  Metric  Conversion  Act  of  1975  and  study 
how  other  countries  made  the  successful  changeover  in  recent 
decades. 
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Figure  3.  Circular  fiber  embedded  in  an  infinite  half-plane  under  an  indentor 


FIBER 

MATRIX -1st  PART 


Figure  4.  Finite  element  model  with  boundary  conditions  and  external 
pressure  acting  on  the  matrix  for  finite  element  results. 
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Figure  5.  Indentor  load  as  a  function  of  axial  relative  displacement  at  fiber- 
matrix  interface  for  flat  indentor  and  uniform  pressure. 
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Normalized  Indentor  Load,  (P/mb^) 
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Figure  6.  Indentor  load  as  a  function  of  axial  relative  displacement  at  fiber- 
matrix  interface  for  uniform  pressure. 
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Normalized  Indentor  Load,  (P/  tc  b^) 
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Figure  7.  Indentor  load  as  a  function  of  axial  relative  displacement  at  fiber- 
matrix  interface  for  transversely  isotropic  fiber. 
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Before  Loading 


At  Peak  Load  After  Unloading 


Figure  8.  Schematic  of  the  slice  compression  test. 
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Axis  of  Symmetry 


Figure  9.  Boundary  conditions  and  geometry  for  analytical  models. 
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Figure  10.  Variation  of  debond  length  with  fiber  volume  fraction  (net 
load=150  MPA,  friction  factor  =  0.01) 
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Maximum  Protrusio 


Fiber  Volume  Fraction 


Figure  11.  Variation  of  maximum  protrusion  with  fiber  volume  fraction  (net 

load=150  MPa,  friction  factor  =  0.01) 
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Residual  Protrusion  (jit 


Figure  12.Variation  of  residual  protrusion  with  fiber  volume  fraction  (net 
load=150  MPa,  friction  factor  =  0.01) 
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Radial  Distance  (|J.m) 


Figure  14.  Stress  distribution  over  fiber  and  matrix  surface  for  a  load  of  150 
MPa  over  the  aluminum  plate  (fiber  radius  =  0.01  mm,  fiber  volume  fraction 

=0.32) 
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